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ABSTRACT 

There is a sign that long-duration gamma-ray bursts (GRBs) originate from 
the core collapse of massive stars. During a jet puncturing through the progenitor 
envelope, high energy neutrinos can be produced by the reverse shock formed at 
the jet head. It is suggested that low-luminosity GRBs (LL-GRBs) are possible 
candidates of this high energy neutrino precursor up to ~ PeV. Before leaving 
the progenitor, these high energy neutrinos must oscillate from one flavor to 
another with matter effect in the envelope. Under the assumption of a power-law 
stellar envelope density prohle p oc with an index a, we study the properties 
of TeV — PeV neutrino oscillation. We hnd that adiabatic conversion is violated 
for these neutrinos so we do certain calibration of level crossing effect. The 
resonance condition is reached for different energies at different radii. We notice 
that the effective mixing angles in matter for PeV neutrinos are close to zero so the 
transition probabilities from one flavor to another are almost invariant for PeV 
neutrinos. We plot all the transition probabilities versus energy of TeV — PeV 
neutrinos from the birth place to the surface of the progenitor. With an initial 
flavor ratio 0°^ : 0°^ : 0°^ = 1 : 2 : 0, we plot how the flavor ratio evolves 
with energy and distance when neutrinos are still in the envelope, and further 
get the ratio when they reach the Earth. For PeV neutrinos, the ratio is always 

Earth. In addition, we discuss the 
dependence of the flavor ratio on energy and a and get a pretty good result. 

This dependence may provide a promising probe of the progenitor structure. 

Subject headings: gamma-ray bursts: general — neutrinos 
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1. Introduction 

The idea that GRBs can serve as the sources of high energy neutrinos has long been dis¬ 
cussed (Waxman & Bahcall 1997, 2000; Dai & Lu 2001; Li et ah 2002; Dernier et ah 2003; 
Razzaque 2013; Vietri 1995; Meszaros & Rees 2000; Guetta et ah 2004; Murase & Nagataki 
2006; Rees & Meszaros 2005; Murase 2008; Wang &: Dai 2009; Gao et ah 2012; Meszaros & Waxman 
2001; Enberg et al. 2009; Murase & loka 2013; Pruet 2003; Razzaque et ah 2003; Ando & Beacom 
2005; Horiuchi & Ando 2008). Before breaking out, a relativistic jet punctures through 
the stellar envelope and transports energy to electrons and protons via shock acceleration 
(Zhang & Meszaros 2004; Piran 2005; Meszaros 2006; Woosley 1993; MacFadyen & Woosley 
1999). Accelerated electrons dominate radiation by synchrotron or inverse Gompton mech¬ 
anism while accelerated protons produce neutrinos by proton-proton collision and photo- 
pion process (Waxman & Bahcall 1997, 1999; Rachen & Meszaros 1998; Alvarez-Muniz et al. 

2000; Bahcall & Waxman 2001; Guetta & Granot 2003; Murase et al. 2006; Becker 2008). 

This neutrino signal is prior to the main burst, as a precursor. However, it has been argued 
that this high energy neutrino precursor of energy ranging from TeV up to PeV cannot be 
produced for a typical GRB (Levinson & Bromberg 2008; Katz et al. 2010; Murase & loka 
2013). The reason is that the reverse shock occurring at the interface of the jet head would 
be radiation mediated, resulting in an inefficient shock acceleration. Nevertheless, a matter 
is different for LL-GRBs (Murase & loka 2013; Xiao & Dai 2014). Due to its low power, the 
Thomson optical depth is low even inside a star so that efficient shock acceleration would be 
expected. We assume the same LL-GRB as our previous work (Xiao & Dai 2014), in which 
we have shown that our LL-GRB is responsible for TeV — PeV neutrinos. Further in this 
paper we focus on the oscillation properties of these high energy neutrinos. 

Neutrino oscillation in matter has been studied for a long while. The resonant con¬ 
version of neutrinos from one flavor to another with matter effect was hrst discussed in 
solar neutrino problem (Ghen 1985). While propagating in a medium, Vf. interacts via neu¬ 
tral current (NG) and charged current (GG), whereas and interacts only via NG. This 
mechanism is called the Mikheyev-Smirnov-Wolfenstein (MSW) effect (Mikheyev & Smirnov 
1985; Wolfenstein 1978). As we know, solar neutrinos are relative low in energy (< lOOMeV), 
as is the same for supernova neutrinos, of which the research is becoming more and more 
attractive (Lund & Kneller 2013; Scholberg 2012). In the meantime, osillation of GRB neu¬ 
trinos is yet less understood. With energy up to PeV, the oscillation properties of GRB 
neutrinos are different compared to MeV neutrinos. Kashti & Waxman (2005) discussed the 
electromagnetic and adiabatic energy losses of vr's and p's which would modify the flavor 
ratio produced by GRBs and concluded that the flavor ratio on Earth is 

1 : 1 : 1 at low energy to 1 : 1.8 : 1.8 at high energy with a transition energy around lOOTeV. 

Mena et al. (2007) discussed high-energy neutrinos produced in optically thick astrophysical 
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objects in the energy range 0.1 — lOOTeV. They studied in detail the shower-to-muon track 
ratio R and discussed its variation with source properties and neutrino oscillation param¬ 
eters. Razzaque & Smirnov (2010) also presented a detailed and comprehensive study of 
flavor conversion of neutrinos from hidden sources (jets) but their results differ from the 
results of Mena et al. (2007). Sahu & Zhang (2010) showed that the resonant oscillation 
could take place within the inner high density region of the choked jet progenitor and the 
hnal flavor ratio detected on Earth is further modihed to either 1 : 1.095 : 1.095 for the 
large mixing angle solution to the solar neutrino data, or 1 : 1.3 : 1.3 for maximal mixing 
among the muon and tan neutrinos in vacuum. Osorio Oliveros et al. (2013) studied the 
three flavor neutrino oscillation of choked GRBs in GeV-TeV energy range for three different 
presupernova star models. They found that for neutrino energies below < lOTeV the flux 
ratio did not amount to 1 : 1 : 1, whereas for higher energy neutrinos it did. Recently, Fraija 
(2014) carried out an analysis of resonance conditions for GeV — PeV neutrinos created in 
internal shocks at different places in the star, estimating the flavor ratios on Earth. 

The main difference of our paper with previous works is at the starting point that we take 
a LL-GRB as the source of high energy neutrinos and hence we focus on higher energy from 
TeV to PeV. We hnd that the mixing angles in matter for PeV neutrinos are close to zero so 
the transition probabilities from one flavor to another are almost invariant, which are different 
with MeV — TeV neutrinos. Thus we get a constant ratio of (f)^^ : (j)^^ : (f)^^ ~ 0.30 : 0.37 : 0.33 
for PeV neutrinos on Earth. Moreover, we discuss the dependence of the flavor ratio on 
the index a and neutrino energy, providing a promising way to probe the GRB progenitor 
structure through a neutrino precursor signal in the future. 

This paper is organised as follows. We present all results in section 2. Subsection 2.1 
is about the density prohle of the envelope and subsection 2.2 is about adiabatic violation 
for neutrinos above TeV. In subsection 2.3 we discuss the three neutrino mixing both in 
the envelope and in vacuum, and then the dependence on envelope density prohle index is 
exhibited in subsection 2.4. We hnish with discussions and conclusions in section 3. 


2. Neutrino Mixing 
2.1. Density Profile of the Envelope 

In this subsection, we take a = 2 as our premise, and discuss the dependence on a later 
in subsection 2.4. 

We assume a power-law envelope density prohle p(r) = Ar~°', where A = (3—a)MHe/(47ri?^““) 
and 2 < a < 3 with Mne and R being the mass and radius of the helium envelope. For a 
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helium core of mass Mne = 2Mq and radius i? = 4 x 10^^ cm, the ambient envelope den¬ 
sity can be expressed as p(r) = 7.96 x gcm“^. The number density of electrons in 

the envelope is iVe(r) = where the number of electrons per nucleon 17 needs to be 

obtained. 


The Saha’s equation reads 


log 


Nr+l 

Nr 


log 


2n^+i(T) 

UriT) 


c; "indn 

+ -logT-^X,-logPe-0.48, 


( 1 ) 


where Nr,Ur, and Xr stand for the number density, partition function and ionization energy 
of Tth ionization ions, respectively. T is the temperature and = N^kT is the electron 
pressure. For our pure helium envelope, we can get 


N{Ee'^+) 

iV(He+) 


1.65 X 10"^^ 


iV(He+) 
’ iV(He) 


0.69, 


( 2 ) 


if we adopt typical values as T = 15000K, xo = 24.58eV, Xi = 54.4eV. We can see that the 
2nd ionization of helium is negligible so that 17 = jv(H^)^+w(He+) ~ 0.408. 

The effective potential of neutrinos can be expressed as 


feff — y/^CpNe, 


( 3 ) 


where Gp is the Fermi coupling constant. 


2.2. Adiabatic Conversion Violation 


Before we start to consider the neutrino oscillation, we need to check whether the 
adiabatic approximation is valid. The adiabatic parameter 7 is defined as 


7 = 


6m‘^ 


sin 20 tan 26* 


1 


irflnVel ’ 
I fir Ires 


( 4 ) 


where 6nP is the mass square difference between the neutrino mass eigenstates, E is the 
neutrino energy and 9 is the mixing angle. The subscript “res” represents the place at 
which resonance happens. We can easily see that 7 is in proportion to 1/E and the adia¬ 
batic approximation requirement 7 3> 1 is usually fulhlled for neutrinos with relative low 
energy of < lOOMeV, such as solar neutrinos and supernovae neutrinos. However, in this 
paper we focus on high energy from TeV to PeV and we find that adiabatic conversion 
is not applicable now. We plot the adiabatic parameter versus neutrino energy in Fig¬ 
ure 1 and it is obvious that 7 3> 1 is violated for high energy neutrinos. Only 712 and 
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7 i 3 are needed because there are at most two level crossing for neutrinos in the three fla¬ 
vor case (Dighe & Smirnov 2000; Yasuda 2014). The vacuum oscillation parameters we 
adopt are = '^•^4 x 10“^eV^, = 2.43 x 10“^eV^, sin^ 6*12 = 0.307, sin^ 6*13 = 

2.41 X 10 sin^ 6*23 = 3.86 x 10 CP violation phase 6 = I.OStt and normal mass hierarchy 
is assumed (Fogli et ah 2012). 


For the reason above, we are obliged to do certain calibration of level crossing effect. 
The jumping probability is approximately computed by the WKB method (Yasuda 2014; 
Kuo 1989): 


P = 


exp[-f 7 F] -exp[-f 7 ;^] 


(5) 


l-exp[-f7-fe] 

where F is a factor depending on the density prohle. Then in Figure 2 we plot Ph and Pl 


versus neutrino energy, representing the jumping probability from energy eigenstate 
and from z/i^ to V 2 m respectively. 


to 


2.3. Three-neutrino Mixing 

2.3.1. Neutrino Oscillation in the Envelope 


The evolution equation for neutrinos in matter is given by 

i— = [UHoU^ + ( 6 ) 

where Hq = ^diag(—0, 5m\2) and = (z/g, z/^, i^r) is the flavor eigenstate (Fraija 
2014). U is the three neutrino mixing matrix. 


U 


C 12 C 13 

-S12C23 - Ci2S23Sl3e*^ 
S12S23 - Ci2C23Sl3e*'^ 


S12C13 SlsC” 

C12C23 — Si2S23Sl3e®'^ S23C13 

—C12S23 — Si2C23Si3e*'^ C23C13 


(7) 


where Sij = sin 6ij, Cij = cos 6ij. 

Neutrino mixing angles in matter can be expressed as (Fraija 2014; Yasuda 2014) 

• oQ sin 2013 

sm 20 i 3 ,,^ - 


sin 2012 ,m = 


y/(cos 20 i 3 

- 2 FKfr/< 5 m? 3)2 +(sin 2013)2' 


sin 2012 


a/(cos 20 i 2 - 2 EVescl^/Smj 2 y + (sin 2012 ) 2 ’ 


( 8 ) 
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The effective mixing angles 6 * 13 ^^, ^ 12 ,m become maximum 7 r /4 as resonance conditions cos 26*13 
2 £'l/eff/hm ^3 and cos 26*12 = Srril^ are fulfilled respectively. We can see that 

^ 13 ,m, 6 * 12 ,m are functions of neutrino energy and radius since Ks = We plot 6*13 6 * 12 ,m 

versus energy at radius tq = 1.0 x lO^^cm in Figure 3(a) and hnd that ^ 13 ,™, reaches max¬ 
imum for ~ 2TeV neutrinos. Here the radius tq is treated as the born site of high en¬ 
ergy neutrinos. The reason is that the Thomson optical depth is tt = naxl = —cft^ — 

0.158( 3Q42tgs-i )^(^)^(^)^( 4xioiicm )~^ at r = ro we have tt ^ 1, thus ensuring 
the efficient shock acceleration (Xiao & Dai 2014). We also plot the effective mixing angles 
versus energy at progenitor surface i? = 4 x lO^^cm in Figure 3(b) and how they evolve with 
propagation distance for ITeV neutrinos in Figure 3(c) respectively. In addition, we kindly 
hnd that effective mixing angles tend to be constant zero for PeV neutrinos, which is the 
main reason for constant havor ratio and will be shown later. 


The transition probability from one havor to another after level crossing calibration can 
be expressed as (Yasuda 2014) 


/1-Pi Pi 0\ 

P{ya^ 1 ^ 0 ) = \Up2,m? Pl 1 - Pl 0 

Vo 0 1/ 


/ 1-Ph 0 Ph \ / \U^i,m\^ \ 

X 0 1 0 |Pa2,m/ 

V Ph 0 I-Ph J \ |Pa3,m/ / 


where a, /9 ^ e, p., r and Uai,mi U^j^m are mixing matrix elements in matter. 


( 9 ) 


We plot the nine mutual transition probabilities between three neutrino havors as func¬ 
tions of neutrino energy at diherent radii in Figure 4(a)- 4(b) and the evolution with 
propagation distance for given energies in Figure 4(c)- 4(d). Given an initial havor ratio 
= 1 : 2 : 0 , we can plot how the havor ratio changes with energy and distance 
when neutrinos are still in the envelope in Figure 5(a)- 5(d). The deviation from 1 : 2 : 0 at 
born site tq is due to level crossing ehect for diherent energies. We can see the trend that 
the havor ratio changes more gently for neutrinos with higher energy. For PeV neutrinos, 
the havor ratio is almost a constant value ~ 0.21 : 0.77 : 0.02 in the envelope. 


2.3.2. Neutrino Oscillation from Progenitor to Earth 

Neutrinos go through vacuum oscillation after leaving the progenitor surface and has 
been well understood. We can express the transition probability P°^ as the hrst-order 
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expansion of the small parameter sin 6^13 (Xing & Zhou 2006): 

= 1 - ^ sin^ 2012 , 

Pe^l = ^ sin^ 2012 COS^ 023 + ^ siu 4012 sin 2023 sin 0i3 cos S , 

= P°g = ^ sin^ 2012 sin^ 623 - ^ sin 40i2 sin 2023 sin 0i3 cos 6 , 

= 1 — ^ sin^ 2023 ~ ^ sin^ 20i2 cos^ 623 — ^ sin 40 i 2 sin 2023 cos^ 623 sin 0i3 cos 6 , 

P°^ = P°^ = ^ sin^ 2023 - ^ sin^ 20i2 sin^ 2023 + ^ sin 40 i 2 sin 4023 sin 0i3 cos 6 , 

P)(_ = 1 “ ^ sin^ 2023 ~ ^ sin^ 20i2 sin"^ 023 + ^ sin 40 i 2 sin 2023 sin^ 023 sin 0i3 cos 6 . ( 10 ) 

The flavor ratio on Earth is 



We plot the flavor ratio versus neutrino energy on just leaving the progenitor and on Earth 
in Figure 6(a)- 6(b). We can see clearly that the flavor ratio varies with energy in the range 
of less than lOOTeV, while the flavor ratio keeps invariant : (f)^^ ~ 0.30 : 0.37 : 0.33 

for PeV neutrinos. 


2.4. Dependence on a 

It is reasonable to argue that the hnal neutrino flavor ratio depends on the density 
prohle of the progenitor envelope. Apparently, with the same assumed envelope mass and 
radius, different values of the power law index lead to different ambient envelope densities, 
thus effective potentials are different. This will have an impact on the resonance conditions, 
effective mixing angles in matter and transition probabilities. In this subsection, we would 
like to investigate how large this impact could be. Here we adopt the same helium progenitor 
but with different power law index a = 2.5, 2.7. Respectively, we can write them as p{r) = 
2.52 X 10^®r“^'^gcm“^ and p{r) = 3.16 x 10^®r“^'^gcm“^ and all calculations have been 
repeated for these two cases. 

We present our result in Figure 7. For simplicity, we only show the flavor ratio at the 
surface of the progenitor (Figure 7(a)) and on Earth ( Figure 7(b)). It is clear that a has 


an impact on the flavor ratio. At a given radius, the resonance conditions are shifted to 
higher energies for larger a. The lower limit of neutrino energy for constant flavor ratio is 
highest for a = 2.7, which is several PeV, compared with subPeV for a = 2.5 and ~ lOOTeV 
for a = 2. The reason is that the effective potential of neutrinos is lower for a steeper 
envelope density profile at the same radius, so higher neutrino energies are required to reach 
resonance conditions. Furthermore, the flavor ratio is evidently different before it reaches 
constant for the three cases, so that we can use the observed ratio of neutrino energy from 
TeV to several hundred TeVs to probe the stellar structure, provided that we can observe 
precursor neutrinos of a GRB with km^ scale detectors like IceCube in the future. 


3. Discussions and Conclusions 

High-energy neutrinos can be produced while the jet is still propagating in the envelope 
and LL-GRBs with typical parameters are responsible for TeV — PeV neutrinos. These 
neutrinos will oscillate with matter effect in the envelope and go through vacuum oscillation 
after leaving the progenitor till they arrive at the Earth. We investigate the three-neutrino 
mixing properties with matter effect and then get an expected flavor ratio on Earth, given 
an initial ratio 0°^ : 0°^ : 0°^ = 1:2:0. 

We notice that adiabatic conversion is violated because level crossing effect is non- 
negligible for such high energy neutrinos. After calibrating this effect, we get the neutrino 
mixing angles in matter and nine transition probabilities. We And that the effective mixing 
angles tend to be zero for neutrinos on the high energy end (~ PeV), resulting in constant 
transition probability and constant flavor ratio. For PeV neutrinos, we always get : 

(pi,^ ~ 0.30 : 0.37 : 0.33 on Earth. 

From our expectations, the final neutrino ratio will depend on the density profile param¬ 
eter a. We take a = 2, 2.5, 2.7 to verify the dependence and the result is clear in Figure 7. 
While the flavor ratio on the high energy end keeps constant, it is evidently different for 
neutrinos with energy from TeV to several hundred TeVs , thus providing a promising way 
to probe the stellar structure in the future. 

In this paper, we only investigate the high-energy neutrino oscillation of one typical LL- 
GRB for the given parameters. Ghanging these parameters may have an impact on the flavor 
ratio-neutrino energy dependence in TeV range: at a given radius, the resonance energy may 
differ and the lower limit of neutrino energy for constant flavor ratio is also different, similar 
to the features shown in Figure 7. However, it does not influence the final flavor ratio of 
PeV neutrinos. That is because the envelope is always too “dense” for PeV neutrinos, the 
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effective mixing angles of which always tend to be constant zero. The case is the same 
when changing envelope temperature T. Though the number density of electrons varies for 
different T, the effective mixing angles for PeV neutrinos are always zero and constant the 
flavor ratio 01,^ : 0^^ : ~ 0.30 : 0.37 : 0.33 on Earth is still in expectation. However, 

this constant value may differ from 0.30 : 0.37 : 0.33 due to the uncertainties of vacuum 
oscillation parameters and neutrino mass hierarchy. 

The IceCube experiment has recently reported the observation of 37 high-energy (> 
30TeV) neutrino events, separated into 28 showers and 9 muon tracks, being consistent with 
an extraterrestrial origin (IceCube Collaboration 2014). To correlate with the observation, 
we take the hypothesis that all the IceCube neutrinos are produced by LL-GRBs before 
jet-breakout. Mena et ah (2014) claimed that : (py^ = 1 : 1 : 1 at Earth is disfavored 

at 92% C.L. with the recently released 3-yr data. On one hand, neutrino flavor ratio in the 
detector may be changed by the Earth matter effect (Varela et ah 2014). On the other hand, 
this does not conflict with our conclusion. We just recommend to do an analysis of observed 
shower-to-track ratio in different energy bins when doing data reduction since we know that 
the flavor ratio depends strongly on neutrino energy and an overall 1:1:1 ratio does not 
make any sense. Nevertheless, we expect constant flavor ratio for PeV neutrinos but we have 
observed only three PeV events now. So our result is to be verihed with a larger dataset of 
TeV — PeV neutrinos of IceCube in the future. If this constant value appears in next few 
decades, when we would have observed tens of PeV neutrino events, we can constrain the 
structure of LL-GRB progenitors and vacuum oscillation parameters by exactly measuring 
this value. If there is no sign of such a constant ratio, the most probable reason is that there 
exist other dominate PeV neutrino sources. The hypothesis that all observed neutrinos are 
produced in the jet propagation process of LL-GRBs may be not so complete because they 
may also origin from other cosmic-ray sources like AGNs or they can be produced in other 
stages of a GRB event such as by internal shocks and external shocks. Especially, we hope 
that one day we could observe the neutrino precursor of a GRB event, this neutrino-GRB 
correlation is crucial for our understanding about the structure of the progenitor envelope 
and the jet propagation dynamics. For a complete comparison with the observation, the 
flavor ratio in TeV range of diffuse neutrino background produced by GRBs needs to be 
done and is beyond the scope of this paper. 
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Fig. 1.— Adiabatic parameters versus neutrino energy of our LL-GRB progenitor. The red 
line represents adiabatic parameter for level crossing between 1st and 3rd energy eigenstates, 
while the blue line is that for level crossing between 1st and 2nd energy eigenstates. 



14 



Fig. 2.— Jumping probabilities versus neutrino energy according to adiabatic parameters. 
The red line represents the jumping probability between 1st and 3rd energy eigenstates, while 
the blue line is that for the jumping probability between 1st and 2nd energy eigenstates. 
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(a) Effective mixing angles versus neutrino energy (b) Effective mixing angles versus neutrino energy 
at born site tq = 1 x lO^^cm at progenitor surface i? = 4 x lO^^cm 



(c) Effective mixing angles evolve versus radius for 
ITeV neutrinos 

Fig. 3.— Effective mixing angles in matter for neutrinos with different energies at different 
radii. In all three subhgures above, the red line represents 6^13 ^ and the blue line is 
All angles are measured in radians. 
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E(eV) 


(a) Transition probabilities versus neutrino energy (b) Transition probabilities versus neutrino energy 
at born site tq = 1 x lO^^cm at progenitor surface i? = 4 x lO^^cm 



(c) Transition probabilities evolve versus radius for (d) Transition probabilities evolve versus radius for 
ITeV neutrinos IPeV neutrinos 

Fig. 4.— Transition probabilities in matter for neutrinos with different energies at different 
radii. Nine different colors represent nine mutual transition probabilities and are listed below: 
Pee red, Pe/r ^ orauge, Per green, P^e blue, P^^ -)■ brown, P^^ purple, P^e -t 

black, Pr/r yellow, Prr cyan. These apply for all four subfigures above. 
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(a) Flavor ratio evolves versus radius for ITeV neu- (b) Flavor ratio evolves versus radius for lOTeV 
trinos neutrinos 



1.0x10 


1.5x10“ 2.0x10“ 3.0x10 

r (cm) 


1.0x10^^ 


1.5x10“ 2.0x10“ 3.0x10^^ 

r (cm) 


(c) Flavor ratio evolves versus radius for lOOTeV (d) Flavor ratio evolves versus radius for IPeV 
neutrinos neutrinos 

Fig. 5.— Flavor ratio in matter evolves versus radius for neutrinos with four different 
energies. In all four subfigures above, the red line represents the fraction of Ve, the green 
line is and the blue line is Vr- 
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(a) Flavor ratio versus energy when neutrinos are (b) Flavor ratio versus energy on Earth 
just leaving the progenitor 

Fig. 6.— Neutrino flavor ratio versus energy before and after long distance vacuum oscilla¬ 
tion. Same as Figure 5, the red line represents the fraction of z/g, the green line is Ufj_ and the 
blue line is 
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(a) Comparison of flavor ratios versus energy when (b) Comparison of flavor ratios versus energy on 
neutrinos are just leaving the progenitors for three Earth for three different progenitor envelope power 
different envelope power law indexes law indexes 

Fig. 7.— The dependence of flavor ratios versus energy before and after long distance vacuum 
oscillation on progenitor envelope power law index. The red lines represent the fractions of 
z/g, the green lines are z/^ and the blue lines are Besides, solid lines are responsible for 

a = 2, dotted lines are for a = 2.5 and dashed lines are for a = 2.7. 




















